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Dynamic experiments with Al-W granular/porous composites revealed qualitatively different behavior with 
respect to shear localization depending on bonding between Al particles. Two-dimensional numerical modeling 
was used to explore the mesomechanics of the large strain dynamic deformation in Al-W granular/porous 
composites and explain the experimentally observed differences in shear localization between composites with 
various mesostructures. Specifically, the bonding between the Al particles, the porosity, the roles of the 
relative particle sizes of Al and W, the arrangements of the W particles, and the material properties of Al 
were investigated using numerical calculations. It was demonstrated in simulations that the bonding between 
the soft Al particles facilitated shear localization as seen in the experiments. Numerical calculations and 
experiments revealed that the mechanism of the shear localization in granular composites is mainly due to 
the local high strain flow of soft Al around the rigid W particles causing localized damage accumulation and 
subsequent growth of the meso/macro shear bands/cracks. The rigid W particles were the major geometrical 
factor determining the initiation and propagation of kinked shear bands in the matrix of soft Al particles, 
leaving some areas free of extensive plastic deformation as observed in experiments and numerical calculations. 



I. INTRODUCTION 

The dynamic behavior of granular/porous reactive ma- 
terials has attracted significant attention due to its pos- 
sible practical applications: reactive structural compo- 
nents, reactive fragments, etc. 1 3 The performance re- 
quirements introduce challenging fundamental problems 
such as combining high strength with ability to undergo 
controlled bulk disintegration producing small sized re- 
active fragments. 

The material structure which can meet this and other 
requirements should be tailored and optimized at the 
mesoscale to produce the desirable mechanical proper- 
ties while still facilitating the release of chemical en- 
ergy. Mesostructural parameters like particle size and 
morphology can affect the strength and shock-sensitivity. 
This can be seen in pressed explosives 4 5 and in reac- 
tive materials like Al-PTFE composites P Additionally, 
mesoscale features like force-chains in granular energetic 
materials may also serve as ignition sites 

Quasi-static, Hopkinson bar, and drop-weight exper- 
iments were performed for PTFE-A1-W and for Al-W 
composites'^ with different particle sizes of W, poros- 
ity, and morphology. W particles were used to increase 
sample density and generate the desirable mode of dis- 
integration of the sample into small sized debris. Multi- 
material Eulerian hydrocode simulations of the dynamic 
tests for the various types of samples were use d to eluci- 
date the observed experimental results.^ 1 16 * 18 1 

The complexity of the dynamic behavior of granu- 
lar/porous materials is caused by the significant role 
of the mesoscale parameters including the particle sizes 



of components, the formation of force chains, the mor- 
phology of particles, and the bonding between parti- 
cles. In these materials, shear localization and fracture 
can be delayed by strain hardening mechanisms such 
as compaction which results in porosity reduction! 17 1 18 1 
At the same time, mesoscale fracture of brittle particles 
can act as a softening mechanism increasing macroshear 
instability! 19 * 20 * In this paper, the dynamic behavior of 
Al-W granular /porous composites, their susceptibility to 
shear localization, and their subsequent fracture at large 
strains was investigated using experiments and numeri- 
cal simulations. Specifically, the bonding between the Al 
particles, the initial porosity, the relative particle sizes 
of Al and W, the arrangements of the W particles, and 
the constitutive behavior of Al was examined in Al-W 
granular /porous composites. 



II. EXPERIMENTS 

High density Al-W granular/porous composites were 
prepared from an elemental powder of Al (Alfa Aesar,- 
325 mesh), and W wires (A-M System, diameter of 200 
/im and length of 4mm) using (a) cold isostatic pressing 
(CIPing) and (b) CIPing followed by vacuum encapsu- 
lation and subsequent hot isostatic pressing (HIPing) to 
create metal bonding between Al particles. All samples 
had the same mass ratios of the Al and W components 
(23.8% Al and 76.2% W, by weight, corresponding to a 
volume ratio of 69.0% Al and 31.0% W) with a theoreti- 
cal density of 7.8 g/cm3. Mixtures of Al powders and W 
(short wires) were placed in a cylindrical stainless steel 
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TABLE I. Initial dimensions of experimental sample 



Sample 


Height [cm] 


Diameter [cm] 


Ciped only 


1.45 


1.58 


Ciped + Hiped 


1.39 


1.54 


Solid Al 6061-T6 


1.25 


1.18 



mandrel with moving pistons and encapsulated in a rub- 
ber jacket providing axial loading during pressurization 
in the CIP chamber. This allowed the preparation of 
samples with the exact sizes and cylindrical shape that 
were necessary for the strength measurements. All sam- 
ples were CIPed at 345 MPa under room temperature 
for 5 min. Subsequent HIPing for some samples were 
carried out at 200 MPa at a temperature of 500 degree 
Celsius with a soaking time 20 min. The density of the 
samples was measured by the hydrostatic method. The 
samples had an average density of 6.8 g/cm 3 after CIPing 
and 7.5 g/cm 3 after CIPing and HIPing. Their dynamic 
behavior and fracture was investigated at a strain rate 
of 1000 1/s under drop weight tests (nominal velocity 
of the drop weight was 10 m/s) using a high velocity 
DYNATUP model 9250HV with an in-house modified 
anvil supporting the sample. The undeformed CIPed 
only, CIPed+HIPed, and solid A16061-T6 samples are 
presented in Fig. [I] After the dynamic tests, shear 
macrocracks were well developed in the CIPed+HIPed 
samples (see Fig. [5] (b)), while no shear macrocracks 
were observed in the CIPed only samples (see Fig. [2] (a)) 
at similar global strains. In the CIPed only sample, Al 
particles were ejected from the outer regions of samples 
leaving areas of interconnected and mechanically locked 
W wires (see Fig. [2] (a) and (e)). The microscopic images 
of the deformed samples representing areas at the vicinity 
of macroscopic shear localization, plastic flow and defor- 
mation of W and Al particles are presented in Fig. [3] It 
is evident that the bonding between the soft Al particles 
facilitated shear localization in CIPed+HIPed samples. 
As a result of localized plastic flow in the matrix of Al 
particles we have observed regions of heavily deformed 
and practically undeformed Al particles. 

For comparison, dynamic tests were performed with 
as is A16061-T6 and annealed A16061-T6 (temperature 
425C, 2.5 hours in vacuum) cylindrical samples under 
similar conditions of impact and similar or larger final 
strains in comparison with CIPed or CIPed+HIPed sam- 
ples (see Fig. [2}. The as is A16061-T6 and annealed 
A16061-T6 samples did not exhibit signs of shear local- 
ization at this level of strains. 



III. NUMERICAL MODELING AND DISCUSSION 

A number of sample properties were examined to see 
how they influenced shear instability and localization 
during a dynamic test. The role of the metallic bonding 




FIG. 1. (Color online) Initial CIPed (a) and CIPed+HIPed 
Al-W samples (b) and a sample of solid 6061-T6 (c). Sample 
diameters and heights are given in Table [I] Mesostructures of 
CIPed and CIPed+HIPed samples before dynamic deforma- 
tion are presented in (d) and (e) respectively. 




(e) 

FIG. 2. (Color online) Deformed shapes of the CIPed (a) and 
CIPed+HIPed samples (b). For comparison, the deformed 
annealed A1-6061-T6 samples and A1-6061-T6 are presented 
correspondingly in (c) and (d). A high speed camera photo 
(e) captures the behavior of Al fragments as they are ejected 
from the CIPed sample during a dynamic test. Arrows are 
added to show areas where the Al particles on the outer area 
of the CIPed sample have been ejected (a) and the shear bands 
that form in the CIPed + HIPed samples (b). 



(a) (b) 

FIG. 3. (Color online) Microstructure of deformed 

CIPed+HIPed samples (a) and (b) illustrating a region of 
largely undeformed Al particles neighboring heavily deformed 
and fractured areas of Al particles with different arrangements 
of less deformed W rods. 



TABLE II. Particle sizes and distributions 





Particle Diameter 


% of Material 




//m 




-325 mesh 


32 


50 




15 


40 




4 


10 


Fine 


40 




Coarse 


200 





was examined in an attempt to understand a difference 
between dynamic behavior of CIPed samples (no metallic 
bonding between Al particles) and CIPed+HIPed sam- 
ples (metallic bonding between Al particles). The roles 
of the initial porosity, the relative particle size of Al and 
W, the initial arrangement of the W particles, the con- 
finement, and the constitutive behavior of the Al were 
also explored to see their effect on the shear instability 
in the sample. A two-dimensional Eulerian hydocode^ 
was used to simulate the behavior of the samples dur- 
ing the drop weight tests. Due to the smallness of the 
particles relative to the global sample size, small repre- 
sentative elements of the microstructure were used in all 
the simulations. 



TABLE III. Material volume fractions and porosity in each 
of the simulated samples 



Sample in 


Volume 


Volume 


% Initial 


Coresponding 


Fraction 


Fraction 


Porosity 


Figure 


Al 


W 




l.a 


0.72 


0.28 


8.3 


l.b 


0.74 


0.26 


0.0 


2 


0.67 


0.33 


8.5 



(a) (b) 

FIG. 4. (Color online) Initial arrangements of the Al (blue) 
and the W (red) particles in the Al-W granular/porous repre- 
sentative samples with similar sized Al an d W particles cor- 
responding to -325 mesh (see Tables |lT| and |III[ ) with .3% (a) 
and 0% (b) initial porosity. Representative element sample 
sizes are 0.02 cm (horizontal) by 0.015 cm (vertical). These 
samples were used for both the bonded and the unbonded 
simulations. 




FIG. 5. (Color online) Initial arrangement of the Al (blue) 
and the W (red) particles in the Al-W granular/porous sample 
with fine Al (40 m) and coarse W (200 m) particles, and a 
porosity of 8.5%. This arrangement was used both for the 
bonded and the unbonded Al particles. Sample size is 0.0625 
cm (horizontal) by 0.1 cm (vertical). 



A. The initial particle arrangements, the constitutive 
model, and the boundary conditions 

The initial particle arrangements of W and Al used in 
the numerical simulations are shown in Fig. [4] through 
Fig(6] The description of the particles sizes for both the 
Al and W are presented in Tables [TT] and |III| . Samples 
in Fig. [4] (a) and (b) were used to investigate the role 
of porosity and metallic bonding on the shear instabil- 
ity and formation of the shear bands in both the bonded 
(CIP+HIP) and unbonded (CIP only) samples with sim- 
ilar sized Al and W particles. The sample in Fig. [5] was 
used to investigate how relative sizes of W and Al affected 
the shear instability in both the bonded and unbonded 
samples. Additionally, the sample in Fig(5]was used to 
examine cases with confinement boundary conditions as 
well as the variations for simulations looking at the role 
of constitutive behavior in Al. Samples in Fig. [6] were 
used to investigate the role of initial mesostructure on 
the shear instability in the bonded case. 

A standard Johnson- Cook with failure^ material 
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Velocity BC 



Velocity BC 



FIG. 6. (Color online) The initial arrangement of the 200 
m diameter W (red) particles embedded in a fully dense (0% 
porosity) Al (blue) matrix. Samples have sizes of 0.3 cm (hor- 
izontal) by 0.3 cm (vertical). Both samples have the same 
volume fractions of Al and W of 0.70 and 0.30 respectively. 



TABLE IV. Consitutive model and equation of state material 
parameters 





Al 6061-Td^ 


y^23l25l 


Johnson-cook 



p [gm cm x ] 


2.7 


16.98 


G [Mbar] 


0.26 


1.24 


A [Mbar] 


2.24-10" 3 


1.506-10" 2 


B [Mbar] 


1.114-10" 3 


1.765-10~ 3 


n 


0.42 


0.12 


c 


2.0-10" 3 


1.60-10" 3 


m 


1.34 


1.0 


C p [J gm" 1 K- 1 ] 


0.89 


0.13 


Tmelt [K] 


930 


1728 


T r0 om [K] 


300 


300 


Di 


-.077 


0.0 


D 2 


1.45 


0.33 


D 3 


-0.47 


-1.5 


D 4 


0.0 


0.0 


D 5 


1.6 


0.0 


Mie-Grneisen 


Co [cm /is -1 ] 


0.52 


0.40 


Sl 


1.4 


1.24 


S2 


0.0 


0.0 


S3 


0.0 


0.0 


7o 


1.97 


1.67 



model and the Mie-Grneisen equation of state used for 
both the Al and the W in the simulations. Failure in 
the material is determined through the equivalent plas- 
tic strain as shown in Ref. |22j The material constants 
used in this paper are presented in Table |IV| 

To simulate the drop weight test, a kinematic bound- 
ary condition with constant downward velocity was im- 
posed on the top boundary corresponding to the impact 
speed of the falling weight in the experiment (10 m/s). 
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FIG. 7. Boundary conditions imposed on (a) corner samples 
and (b) internal confined samples. 



For the cases where the small representative element is 
near the outer edge of the large cylinder, the side wall 
boundary conditions shown in Fig. [7] (a) were used. In 
cases where the small representative element is near the 
middle of the large cylinder, the side wall boundary con- 
ditions shown in Fig. [7] (b) were used to account for the 
confinement. A variety of confinement conditions were 
considered in this paper and will be detailed during the 
discussion of Fig. [l3j With the confinement boundary 
conditions, the main concern is keeping the correct de- 
scription of the global deformation of the sample, but due 
to the artificial conditions, the agreement with global 
stresses is sacrificed. This confined geometry produces 
compressive stresses in particles leading to elevated pres- 
sures in the sample. 



B. The role of bonding and initial porosity in samples 
with similar sized Al and W particles (-325 mesh) 

The results of the numerical simulations correspond- 
ing to the initial mesostructures shown in Fig. [4] are 
presented in Fig. [8] corresponding to the deformed sam- 
ple with a global strain of 0.45. The sample with the 
bonded particles and initial porosity of 8.3% began de- 
veloping shear instabilities at the mesoscale once the in 
situ densification that occurred at the early stages of de- 
formation removed a majority of the porosity from the 
sample (the initial average porosity was 8.3% and the 
porosity when shear instability started was 3.2%). The 
bonded sample with zero initial porosity began develop- 
ing shear instabilities on the mesoscale within the first 
5% global strain. The global shear band began devel- 
oping around 0.25 global strain in the case with 8.3% 
initial porosity and around 0.20 global strain in the sam- 
ple with zero initial porosity. This behavior is similar to 
observations in the PTFE-A1-W mixtures and in the ex- 
periments with the Al-W powder mixture P 1 1 1 2 1 1 6 t SH The 
global shear localized zones and the subsequent cracks 
developed at an approximate 45 degree angle but was 
kinked by the rigid W particles resulting in the initial 
shear band forming in a range of angles from 36 to 50 
degrees (see Fig. [9|. The samples in Fig. 8] (a) and (b) 
have a relatively large percentage of small W particles, 
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providing some geometrical homogenization of the mix- 
ture. As a result, the shear bands do not significantly 
deviate from this 45 degree angle. In both samples, the 
mesoscale shear instabilities continue to grow until one 
begins to dominate (at approximately 0.25 global strain) 
and a global shear band develops, which in turn leads to 
the growth of a macro crack. While this simulation only 
looks at a small representative element along the edge, 
the experimental CIP+HIP samples exhibit the same be- 
havior (see in Fig. [2](b)). 

The corner samples that have unbonded particles, sim- 
ulating the CIPed only material, are presented in Fig. 
[8] (c) and (d). These samples do not develop shear in- 
stabilities like the bonded samples with identical initial 
mesostructures. In the unbonded samples, the Al and W 
particles rearrange themselves during the dynamic defor- 
mation, effectively blocking shear instabilities from form- 
ing. Additionally, due to the free boundary on the right 
hand side, the sample undergoes bulk disintegration. Due 
to the two-dimensionality of the simulation and the lack 
of friction in the simulation, the W particles do not of- 
fer much resistance to the movement of the Al out of the 
right boundary. A similar behavior can be seen in the ex- 
perimental samples where the Al on the edge is ejected 
from the sample, disintegrating into agglomerates of Al 
consisting of 5-30 initial sized Al particles. This ejec- 
tion of the Al particles from the network of the W wires 
adjacent to the free boundary was observed during the 
dynamic loading in experiments using a high speed cam- 
era (see Fig. [2] (e)). 

Initial variations in porosity (0% or 8.35%) in the 
bonded sample altered the orientation of the mesoscale 
shear bands despite the practically identical initial 
mesostructure of the W particles. In the sample with 
porosity, the pores are bulk distributed with sizes of ap- 
proximately 4 microns in diameter. This bulk-distributed 
porosity allows for small movements of the W particles at 
the initial stages of deformation due to pore closure. The 
shear instabilities were observed nucleating during this 
early stage of the deformation. The small alterations in 
the W particle arrangement at this early stage therefore 
caused a large deviation in shear bands that developed 
during later stages of the deformation. The unbonded 
samples with and without porosity both showed the same 
characteristic of particle rearrangement blocking shear 
band formation and bulk disintegration. 




(c) (d) 

FIG. 8. (Color online) Role of initial porosity and bonding 
between particles. Deformed mesostructures of -325 mesh Al 
and W particles with the initial arrangement of Al and W 
particles shown in Fig. [4] (a) and (b). Damage is plotted to 
highlight the shear instability and shear bands in the sample. 
All samples shown at global strain of 0.45 with in itial particle 
arrangement presented in: 

(a) Fig. El (a) with bonded particles (CIPed + HIPed). 

(b) Fig. 13 (b) with bonded particles (CIPed + HIPed). 

(c) Fig. 4 (a) without bonded particles (CIPed only). 

(d) Fig. 4 (b) without bonded particles (CIPed only). 




FIG. 9. (Color online) Effect of rigid W particles on the 
kinking of the shear bands in the Al matrix with the W and 
the Al -325 mesh sized particles with the initial mesostructure 
illustrated in Fig. [4] (a). Black dots show the kink locations 
of the shear bands. The angles of the piecewise shear bands 
are 32 (1), 58 (2), 34 (3), 47 (4), 48 (5), and 45 (6) degrees. 



C. The role of bonding and initial porosity on samples 
with fine (40 micron) Al and coarse (200 micron) W 
particles 

Numerical calculations exploring how differences in 
the relative sizes between the W and Al particles (see 
the initial mesostructure in Fig. |5| effect the develop- 
ment of shear instabilities and the formation of shear 
bands/cracks are shown in Fig 10 The W particles in 



crons, the same diameter as the W rods used in the ex- 
perimental samples shown in Fig. [l](a) and (b). The Al 
particles have a diameter of 40 microns, the size of the 
larger particles in the -325 mesh particle size used in the 
experimental sample. The boundary conditions are the 
same as those described in Fig. [7] (a) with a constant 
downward velocity of 10 m/s on the top boundary. 

In the bonded sample simulating the CIPed+ HIPed 
material (Fig. 10 (a) and (b)), shear instabilities be- 



the numerical simulations have a diameter of 200 mi- 



gin developing after the pores close like the porous sam- 
ple with similar sized Al and W particles shown in Fig. 
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[8] (a). Mesoscale shear bands develop at approximately 
0.25 global strain. The shear band that forms is heavily 
influenced (guided) by the W particles in its path. When 
comparing the shear bands in Fig. [lOjwith those in Fig. 
[8j the W particles that are larger relative to the Al have a 
much greater influence on directing and altering the path 
of the shear bands. The larger W particles require the 
shear instability to circumvent a larger radius of the rel- 
atively rigid W particle to connect with instabilities that 
form in the Al on the other side of the W particle. In Fig. 
10 (a) the shear band can be clearly seen circumnavigat- 
ing the W particle in the upper central region, causing 
the shear band to turn almost 90 degrees from one side 
of the particle to the other. The shear band in the lower 
section of the same figure is directed between two adja- 
cent W particles. Due to the larger relative size of the W 
to the Al particles, the heterogeneous nature of the sam- 
ple is increased in comparison to the previous samples 
in Fig. [5J where the small W particles homogenized the 
mesostructure. This creates areas of mostly undeformed 
Al to appear in regions away from W particles and areas 
of deformed Al to appear in regions near W particles. 
Based on the simulations, the path of the shear bands 
closely follows the W particles in CIPed+Hiped mate- 
rial. This behavior is in agreement with the experiments 
(see Fig. (3|. 

In the unbonded sample simulating the CIPed only ma- 
terial in Fig. [To] (c) and (d), the bulk distributed rear- 
rangement of the Al and W particles effectively block the 
development of shear bands. This rearrangement is simi- 
lar to that of the particles in Fig. [8](c) and (d). Like the 
previous case due to limitations in the two-dimensionality 
of the simulations some of the three-dimensional W fiber 
effects cannot be reproduced. However, the trends seen 
in this two-dimensional simulation can also be seen at the 
edge of the experimental sample where the Al is ejected 
from the edge of the cylinder during the dynamic test as 
shown in Fig. [2] (a) and (e). 



D. The relationship between sample strength and porosity 
in corner samples 

The engineering stress at the upper and lower bound- 
aries versus the global strain for the simulation of bonded 
sample described in Fig. 



10 is shown in Fig. 11 The 



solid red curve and the dashed black curve show the av- 
erage engineering stress on the top and bottom bound- 
aries respectively and are nearly identical except for the 
first 4% strain. Based on this, the sample can be seen 
as undergoing a quasi-static deformation. In addition 
to the stress, the percentage of porosity in the sample 
was plotted to show the relationship between the poros- 
ity and sample strength. It can be seen that stresses 
in the sample increased until a global strain of 0.12 was 
reached. The maximum of stress corresponds to the min- 
imum porosity due to in situ densification during the ini- 
tial stage of deformation. After the densification stage 




FIG. 10. (Color online) Role of relative particle sizes of W 
and Al in bonded (a), (b) and unbounded samples (c), (d) 
at two different global strains of 0.30 (a), (b) and 0.45 (c), 
(d). The initial mesostucture is shown in Fig. [5] Damage is 
plotted to highlight the shear instability and the shear band 
development. 



(corresponding to a range of global strain 0.12-0.29), the 
stress begins to drop as the local shear bands begin form- 
ing and the global shear band develops. Later (0.30 
global strain) macrocracking occurred resulting in a de- 
crease in strength accompanied by a rapid increase of 
porosity. This b ehavior is similar to the behavior ob- 
served in! 11 * 12 * 1 ^ ^ The two-dimensional simulation of 
the unbonded material, due to the material rearrange- 
ment, had near zero engineering stress on the top and 
bottom boundary and was not included. The compara- 
ble strength of bonded and unbounded samples in exper- 
iments is mostly due to interconnected W wires that can- 
not be accounted in these two-dimensional simulations. 



E. The role of initial arrangement of W particles in 
CIP+HIP samples (bonded) with zero initial porosity 

Current processing techniques create samples where 
the W fibers are initially placed in a random fashion. 
Multiple randomized samples were created for numer- 
ical analysis to understand how different experimental 
realizations of this randomized W fiber placement ef- 
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FIG. 11. (Color online) Average engineering stress versus 
the global strain at the top boundary (solid red line) and at 
the bottom boundary (dashed black line) of the corner sample 
corresponding to initial mesostructure presented in Fig. [5]and 
velocity of impact 10 m/s. The porosity in the sample versus 
global strain is also plotted (blue dashed and dotted line). 



fects the formation of shear instabilities during dynamic 
deformation. In two characteristic samples with initial 
mesostructure presented in Fig. |6j the W particles were 
randomly placed in an Al matrix such that the volume 
fractions of W and Al in each sample were identical. The 
samples were then subjected to boundary conditions de- 
scribed in Fig. [7] (a) with a constant 10 m/s downward 
velocity. All samples had zero initial porosity to remove 
mesostructural changes due to void closure as seen in 
simulations described in Fig. [8j Results of the two char- 
acteristic simulations are shown in Fig. 12, demonstrat- 
ing variations in the shear instabilities and shear band 
development while highlighting the similarities that the 
samples share. 



Both samples in Fig. 12 have damage plotted to high- 
light the shear band formation, with red corresponding 
to fully damaged Al. It is obvious from the compari- 
son of the deformed samples in Fig. [l2] (a) and (b) that 
changes in the initial mesostructure (see Fig. 6 (a) and 
(b)) greatly alter the location of the shear band despite 
the same volume content of Al and W and the same size 
of W particles. 



While both samples in Fig. 12 have very different shear 
band locations and modes of fracture, there are a number 
of similar characteristics that the samples exhibit. 

First, both samples develop shear instabilities and 
shear bands at approximately the same global strain of 
0.2. It should be mentioned that in these samples with 
zero initial porosity, the global strains corresponding to 
well developed shear localization are similar to the cor- 
responding global strains in the initially porous samples 
(see Fig. |]and Fig. [To]). 

Second, numerous local shear instabilities form due to 



the localized, high strain flow of softer Al around the 
harder W particles. These local instabilities link with 
other local shear instabilities in close proximity creating 
localized shear bands. These bands join with other local- 
ized shear bands until one global shear band transverses 
the entire sample and becomes the dominant macroshear 
band. This macroshear band has a propensity to form 
at a 45 degree angle, with its path locally altered by W 
particles, and spans the entire sample before subsequent 
global shear bands are able to form. This demonstrates 
that the relatively rigid W particles initiate the shear 
instability in these granular composites, enhancing the 
localized high strain plastic flow of Al around the W 
particles. Local plastic strain in the Al around the W 
particles was 3 to 4 times higher than the plastic strain 
in the surrounding Al thus facilitating localized damage 
and subsequent global shear localization within the Al 
matrix. Separate experiments and numerical simulations 
with cylindrical samples made of as is A16061-T6 and 
annealed A16061-T6 did not reveal shear localization at 
similar or larger global strains (see Fig. [2] (c) and (d)). 

Third, the W particles are largely responsible in dictat- 
ing the path of instabilities. This can be seen clearly in 
the bottom center region of Fig. [l2](a), where the shear 
band develops two branches that circumvent a clump of 
W particles and reconvene on the bottom boundary. Ad- 
ditionally, in the top central region of Fig. [l2] (a), a 
closely packed clump of W particles block shear instabil- 
ities from developing in this region. This may be due to 
a similar effect as seen in Fig. [8] (a) and (b) in addition 



to Fig. [101 (a), where having the W particles initially 
aggregated into a close pack causes the shear instability 
to move around the relatively rigid mass into the softer 
area of the surrounding composite. The same effect can 



be seen in all areas of the samples in Fig. 12 where the 
W particles group in very close proximity to each other. 

Finally, due to the random nature of the sample, some 
of the W particles form angular chains with channels of 
Al between them facilitating the growth of shear instabil- 
ities and shear bands along the channel. The sample in 
Fig. [l2] (b) exemplifies this with W particles creating a 
channel structure spanning the entire sample going diago- 
nally from the lower left corner to the upper right corner. 
This channel created a favorable path for the local insta- 
bilities to follow and eventually lead to the formation of 
the shear band in the sample. W particles flow along 
the shear band to help form these angular chains. Rear- 
rangement of the rigid particles inside high strain shear 
flow was observed in explosive compaction on interfaces 
between particles with high strain shear deformation.^ 
In Fig. [l2] (a) a distinct global chain structure of W 
particles is not present like in Fig. [l2] (b), but rather a 
series of shorter chain structures in the left and the up- 
per right region of the sample. This lack of long chains 
may be the reason that the sample accumulated more 
bulk-distributed damage in comparison to Fig. [121(b). 
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FIG. 12. (Color online) Role of initial mesostructure of the 
W particles on the shear localization and the subsequent frag- 
mentation with zero initial porosity. The deformed samples 
correspond to the initial W particles arrangements shown in 
Fig.©: 

(a) at 0.5 strain corresponding to Fig. pi (a) 

(b) at 0.5 strain corresponding to Fig. j6] (b) 



F. The role of kinematic confinement conditions 

Previous simulations examined small representative el- 
ements corresponding to a material on, or very near to, 
the outer edge of the experimental sample. A small rep- 
resentative element near the central region of the sample 
interacts with surrounding material which resists move- 
ment of material leaving the representative element. To 
account for this confinement from the surrounding mate- 
rial, kinematic boundary conditions using a ramped nor- 
mal velocity were imposed on the side. Multiple velocity 
ramps were explored in an attempt to model the large 
variance in the local inhomogeneity within the sample, 
resulting in various local confinements. The first ramp 
velocity tested kept the global Poisson ratio of the small 
representative sample geometrically consistent with that 
of the experimental sample. Two additional constrained 
boundary conditions were tested corresponding to 70% 
and 100% increased horizontal expansion in compari- 
son with the geometrically constrained case. The ver- 
tical strain rate was kept identical for all samples. Fig. 
13| shows the simulation results at 0.30 and 0.50 global 
strains respectively for the samples with the constrained 
side boundary conditions. Damage is plotted to highlight 
the shear instabilities. 

Shear bands developed in all samples at approximately 
0.25 global strain. However, the location and number 
of shear bands differed in each case. Despite these dif- 
ferences, all the samples share a similar characteristic: 
areas of heavily deformed Al particles near the W parti- 
cles, resulting in local damage accumulation and subse- 
quent shear band formation, while material further away 
from the W particles is relatively undeformed. This trend 
seems to be independent of the imposed confinement con- 
ditions. Also, samples in Fig. [l3] (b) and (c) exhibit 
cracks that form between the W particles. The crack lo- 
cations are very similar to the location of cracks in the 
large experimental sample (see Fig. [3](b)). The cracks 
span between the W particles in the areas of the heavily 
deformed Al. In the simulations, the W particles move, 
squeezing the Al particles between the nearby W par- 
ticles, acting as anvils that facilitate the heterogeneous 
deformation of the Al particles. A similar effect may oc- 
cur in the large experimental sample where the Al flows 




FIG. 13. (Color online) Role of different conditions of con- 
finement on shear localization and subsequent fracture. Three 
deformed mesostructures at 0.3 and 0.5 global strains with 
various confinement conditions are presented. All three sam- 
ples have boundary conditions described in Fig. [7] (b) with 
a constant 10 m/s downward velocity imposed on the top 
boundary. Each sample differs only in the imposed horizontal 
velocity on the side boundaries (linear ramp as vertical strain 
goes from to 0.5) in the following way: 

(a) ,(d) A linear ramp from to 7.20 m/s. Corresponding to 
a final horizontal strain equal to 0.58. This produces the rep- 
resentative element to keep the same geometric proportions 
as in the experiments. Sample shown at 0.3 (a) and 0.5 (d) 
vertical strain. 

(b) ,(e) A linear ramp from to 12.24 m/s. This corresponds 
to a horizontal strain 1.0. Sample shown at 0.3 (b) and 0.5 
(e) vertical strain. 

(c) ,(f) A linear ramp from to 14.40 m/s. This corresponds 
to a horizontal strain of 1.16. Sample shown at 0.3 (c) and 
0.5 (f) vertical strain. 



around the W fibers due to the three-dimensional cage 
like network that the packed network of W fibers cre- 
ate. A similar effect of rigid and heavy W particles is 
instrumental also in the shock loading of corresponding 
mixtures allowing tailore d red istribution of internal en- 
ergy between components! 11 * 12 ! 



G. The role of the constitutive behavior of Al 

The constitutive behavior of the Al particles were mod- 
ified to explore how they influenced the shear instabil- 
ity and the shear band formation in the bonded samples 
(CIPed plus HIPed samples in experiments). Simulations 
were carried out with the initial mesostructure shown in 
Fig. [5] and the boundary conditions in Fig. 7 (a) with 
a constant 10 m/s downward velocity on the top bound- 
ary. The initial yield stress of Al was reduced to a very 
low level 20 MPa and the results are shown in Fig. 14 
The sensitivity of shear band development to other pa 
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FIG. 14. (Color online) Effect of the initial yield strength 
on the shear instability in porous granular composites. The 
initial yield strength of the Al is reduced to that of a softer 
Al (20 MPa). The sample is presented with global strains of 
0.30 (a) and 0.40 (b). The damage is plotted to highlight the 
shear bands. 



rameters in the Johnson-Cook model was also explored, 
however, these simulations showed that there were only 
slight changes in the material response due to these al- 
terations and results are not presented in this paper. 

The reduction of the initial yield stress in the mate- 
rial from that of Al 6061-T6 (324 MPa) to that of a 
much softer Al (20 MPa) caused a large change in sam- 
ple response (compare to the results in Fig. 10 (a) and 
(b)). The Al particles in Fig. 



14 



experiences greater 
plastic flow and deformation in all areas of the sample, 
especially in the areas around the W particles. Due to 
the softening of the Al matrix, the development of the 
global shear band was delayed, allowing for more bulk- 
distributed plastic flow of the Al. The reduction of initial 
yield stress in the Al resulted in a shear band develop- 
ing in a different location. This is likely caused by the 
greater movement of W particles in the softer Al altering 
the arrangement of the W altering the mesostructure at 
the earlier stages of sample deformation. 



IV. CONCLUSIONS 

Dynamic experiments with Al-W granular/porous 
composites and numerical simulations revealed the char- 
acteristics that had the greatest effects on shear insta- 
bility and shear band formation. It was shown that in 
simulated CIPed only (unbonded) samples the Al and 
W particles rearranged themselves during the dynamic 
deformation to effectively block shear localization. This 
resulted in the subsequent bulk disintegration of the sam- 
ple in agreement with areas near the outer surface of the 
samples used in the experiments. 

All CIPed+HIPed (bonded) samples exhibited shear 
localization and shear band formation. The shear bands 
nucleated during the initial stages of the deformation in 
Al surrounding the W particles and spread to the nearby 



W particles at angles close to 45 degrees. The shear 
band is kinked by W particles causing the shear band 
path to deviate from the ideal 45 degree angle path dic- 
tated by global geometry. In simulations with relatively 
larger W particles, the path of the shear band was in- 
fluenced to a greater degree than simulations with the 
similar sized W and Al particles due to increased hetero- 
geneity of the sample. It was also shown that variations 
in the initial arrangements of the W particles were the 
main drivers determining where the global shear bands 
formed in the sample. Numerical calculations and ex- 
periments revealed that the mechanism of shear local- 
ization in granular composites is due to a localized high 
strain flow of Al around rigid W particles, causing local 
damage accumulation and a subsequent growth of the 
meso/macro shear bands/cracks. 

A variety of constraint boundary conditions were ex- 
amined to represent the heterogeneous nature of the in- 
ternal structure in the global sample. Each simulation 
showed shear localization occurred between the nearby 
W particles while leaving areas away from the W parti- 
cles relatively undeformed. This result is supported by 
the microstructural features observed in the experimental 
sample. 

Finally, the role of the constitutive behavior in the 
Al was examined in numerical calculations. The results 
showed that a significant reduction of the initial yield 
stress from Al 6061-T6 (324 MPa) to a softer Al (20 
MPa) increased the amount of bulk distributed damage 
and plastic strain in the sample in addition to altering 
the shear band location. 
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